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Abstract: We investigate the exclusive jet mass distribution in e + e~ events, denned with 
a veto on the out-of-jet radiation, at two-loop order. In particular, we calculate the two-loop 
soft function, which is required to describe this distribution in the threshold region. When 
combined with other ingredients using soft-collinear effective theory, this generates the 
complete singular distribution for jet thrust, the sum of the jet masses, at two-loop order. 
The result is in excellent agreement with full QCD. The integrated jet thrust distribution is 
found to depend in an intricate way on both the finite jet cone size, R, and the jet veto scale. 
The result clarifies the structure of the potentially large logarithms (both global and non- 
global) which arise in jet observables for the first time at this order. Somewhat surprisingly, 
we find that, in the small R limit, there is a precise and simple correspondence between 
the non-global contribution to the integrated jet thrust distribution and the previously 
calculated non-global contribution to the integrated hemisphere soft function, including 
subleading terms. This suggests that the small R limit may provide a useful expansion for 
studying other exclusive jet substructure observables. 



1. Introduction 



High energy scattering processes in Quantum Chromodynamics (QCD) often involve the 
production of jets. A QCD observable of particular interest is the jet mass. Jet mass 
plays an important role in many applications, such as in the identification of a boosted W 
boson |l|, |2[ or top quark [||, ||. Jet mass is also important because many event shapes 
studied at e + e~ colliders reduce to some measure of jet mass in particular limits. For 
example, thrust in the dijet limit reduces to the sum of the masses of two hemisphere jets, 
and heavy jet mass reduces to the larger of the two hemisphere jet masses. Thrust and 
heavy jet mass are two of the most well-studied e + e~ event shapes and both have been 
measured with extraordinary precision jS?]. On the theoretical side, next-to-next-to-leading 
order (NNLO) QCD calculations of these observables have been carried out for jets in e + e~ 
events || [7|, [8], ||. Unfortunately, jet physics at hadron colliders is much more difficult, 
both experimentally and theoretically. It is therefore still worth studying e + e~ jets, where 
the analytical calculations can be done, in the hope that universality properties will be 
found which may ultimately make the hadron collider calculations tractable, at least in 
some regime. 

In the limit of small jet mass, one cannot straightforwardly make accurate predic- 
tions for the thrust or heavy jet mass using QCD perturbation theory. In this limit, 
the cross section is dominated by the soft and collinear interactions of massless partons 
which produce large logarithmic contributions of the form a* In-' (m/Q) at each order in 
perturbation theory. Without a correct resummation of these terms to all orders in a s , 
perturbative predictions are not reliable when m/Q —> 0. In recent years, much progress 
has been made on this front and NNNLL resummations of both thrust and heavy jet mass 



are now available [10, 11, O]. These resummations led to a measurement of the strong 



coupling constant, a s , which is competitive with the world average [13]. Soft-collinear 
effective theory (SCET) is the theoretical tool that facilitated many of these important 
developments O, ffR [0], [lTj. Indeed, before the simplifications brought about by SCET, 
thrust and heavy jet mass had only been resummed at the NLL level [18]. 

Although these studies have shown the power of resummation, the results for thrust 
and heavy jet mass are not directly applicable to the study of observables at hadron col- 
liders. At hadron colliders, there is soft and collinear activity in the initial state that 
complicates the definition of jet mass |ll|. In reality, experiments at hadrons colliders 
use inclusive jet definitions, and an explicit veto procedure to reduce the huge QCD back- 
ground. Therefore, it is interesting and non-trivial to try and carry out the resummation 
of a QCD observable in the presence of a jet veto. In fact, some work in this direction has 
already been done and it is known that, in this case, a new class of so-called non-global 
logarithms (NGLs), qualitatively different from the class of logarithms encountered in the 
resummation of thrust and heavy jet mass, appear which also need to be resummed in 
certain kinematical limits [^, pjjfl . Unfortunately, the resummation of these loga- 



rithms is not straightforward (see Refs [|2C|, |22|, 24, 25, |26|]). Over the last couple of years, 



there has been renewed theoretical interest in exclusive jet mass distributions and NGLs, 
both with and without SCET |7|, ||, ||, ||, |l], ||, |3|. 



- 1 - 



In a previous publication, Ref. |3C], a simple observable called jet thrust was intro- 



duced. It was subsequently studied in Refs. [33, 34]. Jet thrust, denoted r w , is defined as 
the sum of the masses of the two hardest jets in an event with a veto that restricts the total 
out-of-jet energy to be less than uj. This observable, being a single variable (once the jet size 
R and the veto scale u are fixed) is easy to study qualitatively with one-dimensional plots. 
It is also free of a subset of non-global logarithms of the form In (7714/771,2). In Ref. |30|, it 
was observed that the soft function for jet thrust seems to refactorize at small R. Heuristic 
arguments for this refactorization were given, as was a convincing numerical comparison of 
predictions derived from this refactorization to the singular part of the exact 0{a^) results 
in QCD derived using the program event2 pq|. The leading non-global logarithms 



for jet thrust, which were ignored in Ref. [p0 |, were subsequently studied in Refs. |33|, 



Understanding analytically how the refactorization conjectured in Ref. [30] works, as well 
as the structure of the subleading non-global contributions, requires the two-loop result. 
The two-loop result is the subject of the present paper. 

In Ref. p4| , the exact O (oq) hemisphere soft function was calculated. This facilitated 
an analytical study of the hemisphere mass distribution at NLO. In particular, subleading 
non-global logarithms were uncovered for the first time, as well as an intricate non-global 
functional dependence, in sharp contrast to the simple powers of logarithms suggested 



by |Tjfl • Parts of this result were independently confirmed by two other groups [31]) ^3 . 

In this paper, we go one step further and perform the calculation of the soft function 
for jet mass with a jet veto at 0{a 2 s ) using jets of cone size R. This produces an exact 
form for the singular part of the differential jet thrust distribution. In fact, we present 
the complete scale-dependent contribution to the integrated jet thrust distribution in the 
threshold region. Our result is exact up to the addition of a t u - and ^-independent function 
of R. This function is calculable, but does not affect the differential distribution, and is 
therefore not essential for our present purposes. Our results allow us to both test the 
refactorization ansatz of Ref. [30] and determine the precise form of the non- global terms 



which arise in the integrated jet thrust distribution. 

A number of authors have already observed that there are similarities between the 
leading NGLs in jet observable distributions defined using narrow anti-A>r jets and the 
NGLs that appear in the hemisphere mass distribution [29, 30, [n], |33|| . The i?-dependent 
coefficient of the leading non-global logarithm 

/Wm'(^) (1.1) 

has been calculated for various algorithms, for arbitrary R. However, all existing NGL 
studies with R dependence focus on the extraction of the leading non-global logarithm, 
which can be done without a complete two-loop calculation. 

In this work, we calculate the complete non-global contribution to the integrated jet 
thrust distribution, up to a scale-independent function of R. Remarkably, in the small 
R limit, we find that this non- global contribution to the integrated t u distribution can 
be expressed in terms of the non-global contribution to the integrated hemisphere soft 
function calculated in Ref. p3. The correspondence actually reproduces the entire non- 
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global functional form, not just the non-global large logarithms. We also confirm that all 
of the singular terms in R in the integrated jet thrust distribution can be reproduced by 
the refactorization formula of Ref. |3(|, provided that one takes into account the natural 
R dependence in the argument of the non-global function. 

This paper is organized as follows. In Section ^, we briefly review the factorization 
theorem for the distribution and recall the refactorization ansatz for the t u soft function 
proposed in Ref. [3C]. In Section 0, we summarize the calculation of the r w soft function, 



at times focusing only on those moments of the integrals which ultimately contribute to 
the scale-dependent terms in the integrated t u distribution. In Section |], we compare the 
prediction of SCET to full QCD and use the difference of our results and the prediction of 
the integrated refactorization ansatz to define the (^-independent) non-global contribution 
to the integrated r w distribution. In Section ||, we take the small R limit of the non-global 
contribution to the integrated t u distribution. We then make precise the correspondence 
alluded to above between this non-global contribution and the non-global contribution to 
the integrated hemisphere soft function. Finally, we analyze the NGLs in detail, discuss 
their R dependence, and present a refined refactorization formula for the r u soft function 
consistent with our two- loop results. In Section ||, we present our conclusions and discuss 
some interesting open problems. In Appendix |A| we present some of the technical details 
of our calculation of the integrated t w distribution and in Appendix [B| we collect two 
rather lengthy analytical expressions that would have been awkward to define in the text 
where they first appear. A Mathematica notebook with these expressions is available upon 
request. 

2. Jet Thrust in SCET 

Soft functions are squared matrix elements of Wilson lines integrated against some mea- 
surement operator. They are integral to QCD resummation studies and have appeared in 



the literature in many different contexts |3(j, ^0], 41, |^]. In Ref. |3(|, an inclusive jet mass 
observable called jet thrust, r w , was discussed, and after making a conjecture for part of 
the two-loop soft function, the 0{a 2 s ) prediction from SCET for the distribution was 
compared with the output of event2. In this section, we briefly review the definition of 
t u and recapitulate the main results of Ref. [Q which sets the stage for the exact results 
calculated in the present paper. 

Jet thrust r w is defined as follows. First, for some multi-jet event at an e + e~ collider, 
cluster the particles into jets using some jet algorithm with size parameter R. Define A 
to be the energy of the radiation not going into the two hardest jets. 1 Then jet thrust is 
defined as the sum of the squared masses of the two hardest jets normalized to the center 
of mass energy if A < uj and as otherwise. In symbols, 

T „ = M±M e( .-A), (2.1, 



1 \ can also be defined as the energy of the third hardest jet, as it was in |30|. The definitions are 
equivalent definition as far as the differential t w distribution at two loops is concerned. 
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where M\ and M2 are the masses of the two hardest jets and 0(x) is the unit step function. 



In the limit t w <C 1 and ui/Q <C 1, the jet thrust distribution factorizes |30 
1 dcr 



00 dr^ 



(Q 2 , fj) jdk L dk R dM 2 L J(Ml - Qk L ,n)J(M 2 R - Qk R , fj,) 

M 2 + Ml 

up to power corrections in r w and uj/Q. In Eq. (2.2), the hard function, H(Q 2 ,fi), and the 
jet function, J(p 2 ,fi), are the same as the ones which appear in the factorization theorems 



d\S R (k L ,k R ,X,^d(r u - L ^ 2 R ) , (2.2) 



for the thrust and heavy jet mass distributions [28, |4q|. In principle, one could use R 



R ij=o l-cos% (2.3) 



dependent jet functions [27, 28, 44|. However, the R dependence of such jet functions 
shows up only in terms non-singular in the jet masses. Therefore, due to the fact that we 
are working in the threshold region where SCET is valid, we are free to use the simpler 
inclusive jet functions. 

The distribution of r w depends on the precise jet definition used. In the QCD calcu- 
lation, the jet definition affects the whole distribution. In the SCET calculation, it only 
affects the soft function. In Ref. [ 30 1 , the Cambridge/Aachen (C/A) algorithm was used. 
The C/A algorithm first calculates the distances between all particles i and j 

1 

2 

and then merges the two closest particles into a single particle by adding their four- 
momenta. The algorithm stops when no two objects are closer than a given jet size, 
R, to one another. The Cambridge/Aachen algorithm has many appealing qualities and is 
infrared-safe. The corresponding t w soft function, however, appears challenging to calculate 
at two loops. 

Fortunately, there is a simple cone algorithm one can use instead. We define cone 
algorithm jet thrust as follows. First, find the thrust axis in the event. Then define the 
distance Ri n (Rm) between particle i in the event and the thrust axis n (n) as 

Rin = ~(1 -COsft n ) Rin = ~(1 -COS ft a), (2.4) 

where ft n (ft a) is the angle between particle i and n (n). If Ri n (Ri n ) is less than R, then 
cluster particle i into the right (left) jet. Because we use the thrust axis, this algorithm is 
infrared-safe at e + e~ colliders. 

The jet thrust factorization formula bears a close resemblance to the factorization 
formula for thrust [10, 17]. The only difference is in the soft function. In the case of thrust, 
the soft function is just the hemisphere one, Shemi(fc.L) &R> A*)> and depends on two scales, 
ki and k R . The t w soft function depends on the scales A and R as well and it is therefore 
significantly more complicated. The soft function is defined as 



S R (k L ,k R ,X,fj,) = ±^(0\Y fl Y n \X a )(X a \YiY*\0) 

x5(k R - n ■ P?)5(k L -n.pL)6(\-E Xs ), (2.5) 
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where (Px) * s ^ ne four- momentum of the soft radiation clustered into the right (left) 
jet, and Ex a is the total energy of the out-of-jet radiation. At zeroth order, the soft function 
is simply a delta function 

S R (k L , k R , A, n) = 6(k L )6(k R )S(X), (2.6) 
reflecting the fact that there is no soft radiation. The first non-trivial corrections come at 



0{ct s ) (what is traditionally called leading-order (LO)) and were calculated in Refs. [28| 
and @. 

Using the renormalization-group (RG) invariance of the factorization formula, one can 
show that the soft function must have the form 

S R {k L , k R , A, fi) = S fl (k L ,fj,)S f _ l (k R , (j,) <8) Sf(X, k L , k R ) , (2.7) 

where <S> denotes a convolution (it is a product in Laplace space). Here S^k,^) is the soft 
evolution kernel that precisely cancels the RG evolution of the jet and hard functions |L2"|, 



17, 37 1 . It follows that the soft function can be written as 

S R {k L ,k R ,X,ti) =S%(k L ,v)S%(k R ,riS R ut (\,fx)®S f R (\,k L ,k R ). (2.8) 

There is not actually much content in this separation until the objects involved are given 
precise definitions. For example, we can take S R nt (X,/i) = 1, thus reducing Eq. ( gj ) to 
Eq. fl2.7j), A refactorization formula like this has appeared multiple times in the litera- 



ture |27], |2^, |3(| 33 1 . However, without a precise statement about S R nt (X,fi), Eq. (^ 
actually has less content than Eq. fl2.7|) . 

In Ref. [3C] it was argued that, in the small R limit, S R (k,/i) and S R ut (X,/i) should 



act like soft functions in their own right, with anomalous dimensions. Only the function 
S R (k,n) can have a /U-dependent anomalous dimension, since the soft-collinear region of 
phase-space lies within the jet. Thus, the anomalous dimension of S R (k,/i) has a cusp 
part which agrees with that of S^k,^). In Ref. [3C] it was conjectured that the regular 



anomalous dimensions of the in and out soft functions, at small R, should be 

1% = \ {is + r cusp In ) and = " r -P ln ( > ( 2 - 9 ) 

where r cusp is the usual cusp anomalous dimension and 75 is the regular anomalous di- 



mension of the thrust soft function [111. Note that this ansatz has the cusp anomalous 
dimension generating something other than m(/i) in the anomalous dimension. The con- 
nection of the cusp anomalous dimension to the i?-dependence at LO is also implicit in 
Refs. [27, 28|. Having S R (kL,fJ.), S R (k R ,fi), and S R nt (X,fi) together with their anomalous 
dimensions allows us define S R (X, ki, k R ) precisely: S R (X, kL, k R ) is everything that is left 
over from refactorization with the anomalous dimensions in Eqs. (|2.9|). 

Using the above ingredients, and setting S R (X,kL,k R ) = 1, Ref. [Q predicted the 
differential r w distribution to 0{ot 2 s ) and explicitly compared the result to the output of 
event2. The predictions made from refactorization come very close to reproducing the 



output of event2 in the regime where the factorization formula of Eq. (2.2) is valid. In 
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Figure 1: The difference between coefficient of (f^) 2 in da/d\a{T u ) in full QCD and in SCET. 
We show on the left the results obtained using the refactorization and T cuap ansatz of Ref. [30 . 
The results shown on the right take into account the precise 0(a 2 s ) results derived in this paper. 
The top row uses the Cambridge/ Aachen algorithm for the QCD calculation, whereas the bottom 
row uses a cone algorithm in QCD, in which the jet axis is taken to be the thrust direction. The 
top- left plot is identical to Figure 6 of Ref. pfj. All plots have R = 0.1 and w = 10 _4 Q. 



the top-left panel of Figure ||, we reproduce Figure 6 from that paper. Figure [T] shows 
the differences between the Cambridge/ Aachen algorithm and the cone algorithm (top and 
bottom rows) and the differences between the prediction of refactorization (left) and the 
precise 0{a 2 s ) results derived in this paper (right). The plots displayed in Figure [l] take 
R = 0.1, uj = 10~ 4 Q, and are separated by color structure. 

The plots all show the difference between the differential r w distribution in full NLO 
QCD and the differential r w distribution in various NLO SCET-based approximations. 
If the singular terms in the QCD distribution are correctly reproduced by the various 
approximations, the plots in Figure [l] should tend to zero as r w — > 0. In all four panels this 
vanishing is qualitatively visible. However, the fact that the bottom-right panel vanishes 
best is a sign that both the jet algorithm and non-global logarithms are important (but 
subleading) effects. For example, the leading non-global logarithm, which was set to zero 
in Ref. should contribute a term of the schematic form ln(-r w ) /r w to the differential 
distribution. Two papers, Refs. p3[ and |34}|, have pointed out importance of these NGLs. 



Such terms, however, were beyond the scope of the analyses carried out in Ref. [ 30 1 . 

To gain a better understanding of the refactorization conjecture and the NGLs, we 
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explicitly calculate the soft function to 0(a 2 ) in the next section, with the understanding 
that our ultimate goal will be to study the scale-dependent contributions to the integrated 
t u distribution. As was observed in Ref. |3lJ, the advantage of approaching the calculation 
in this way is that both the //-dependent terms and NGLs of the soft contributions are 
completely transparent at the level of the integrated distribution. There are three color 
structures at 0{a 2 s ): Cp, CfCa, and C^n/Tp. The C F terms are uniquely determined by 
the Abelian exponentiation theorem |45|, [46| and we therefore only present results for the 
CfCa and CpUfTp color structures. 



3. Calculation of the r„, Soft Function and Moments Thereof 



The soft-function is usually expanded in a s as 

oo n 

S R (k L ,k R ,\,fi) = ^2 (|^) S^\k L ,k R ,X,iJ,). 



n=0 



(3.1) 



(2) 

In this section we summarize the calculation of the two-loop r w soft function, S R , at times 
focusing only on those moments of the soft function which will be needed later on when we 
discuss the jet thrust distribution and its non-global structure. We begin with a review of 
the one-loop soft function, establishing our notation, and then describe the NLO results. 
For readers who are only interested in the interpretation of the results, this section can be 
skipped. 

3.1 LO (one-loop) soft-function 

The tree-level result is trivially given as a product of delta functions, 

S { °\k L ,k R ,\,v) = 8(k L )5(k R )5(\). (3.2) 
At first order in a s , the r w soft function in d = 4 — 2e dimensions is given by the integral 2 



t _ x .a _ ^^2^_ ( V 2 ^ E Y f d d q AiFt{k L ,k R ,X) 



S R 1 >(k L ,k R ,\^) = (47ryC F 



47T 



(2vr) c 



q^q- 



(3.3) 



In Eq. ( |3.3D , q + = n ■ q, q = n ■ q, and F\{kL, k R , A) is the function that implements the 
one-particle phase-space cuts 0: 



F 1 (k L ,k R ,X) = -27ri5(q 2 )e(q°) 



®l5{k R - q+)5(k L )5(X) + e q J(k L - q-)5(k R )S(X) 



+®Lt.S [ X 



5{k L )5{k 



R 



(3.4) 



2 In this paper we suppress the ie prescription. This is appropriate since we are dealing only with real 
radiation. 
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In Eq. 0, 



92 = e 



9| = e{rq + -q-) 



w out 



where 



9 (g- - rq+) 6 (g + - rg~) , 
R 



(3.5) 
(3.6) 



1 — i? 

This last abbreviation occurs frequently throughout the rest of this paper. 

Carrying out the integral in Eq. ([3.3]), we find that, up to and including terms of 
0(a„), the t u , soft function can be written as 



s R (k L ,k R ,x,n) = s%(k L ,n)s%(k R ,iJ,)ST(*,») 



where 



and 



l+r 



1 

l+r 



' S 'R Ut (^>/ x ) - <K A ) + ((l- e ) e r(l-e e )(l+0(2A) 1 + 2E j ( (l+ r ) 2 

e 2 F! (l- £,1 + 6,2-6;^) + (e - 1) (1 + r) 2^1 (-6,6,1-e;^ 
ersFx (l - e, 1 + e, 2 - e; ^) + (e - 1)(1 + r) 2^ (-e, e, 1 - e; ^) | 



(3.7) 
(3.8) 

(3.9) 



are the relevant expressions calculated through 0(a s ) and to all orders in e. In Eq. (|3~9|), 
2-Pi(a, b, c; z) is the Gauss hypergeometric function. Expanding the above expressions in e 
reproduces the results published in Ref. |3(| . 

The form of Eq. fl3.9|) suggests that the dependence of the out-of-jet soft function 
is not simple, even at (D(a s ). On the other hand, we note that the 0{a s ) R dependence 
of Sft(k,fi) is just r e ; it factorizes from the rest of the 0{a s ) expression. It turns out that 



the full answer is simply r € times the result obtained several years ago [10, |TJ] for the LO 
hemisphere soft function. As will be shown below, a similar phenomenon occurs at NLO 
for the part of the in-jet soft function that dominates in the small R limit. Stated more 
precisely, at (D(a™) the contribution with n soft partons clustered into the same jet is equal 
to the analogous contribution to the hemisphere soft function multiplied by a factor of r nt . 
This simple relation has its origin in the fact that the integral representation of this part 
of the 0(a™) r w soft function has nice transformation properties under the rescaling 



n. 



Vrn'u 



n. 



-n, 



(3.10) 



This observation will be discussed in more detail in Section |E[ 
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3.2 NLO (two-loop) soft function 

In this section we present the calculation of the NLO t w soft function, focusing at times on 
those moments of the soft function which we will need in subsequent sections. A moment's 
thought reveals that the squared matrix elements written down for the hemisphere soft 
function in Ref. [24] can be used for the r w soft function as well. The two-particle phase- 
space cuts, however, are significantly more complicated. As in the hemisphere case, it is 
convenient to split the calculation up according to how many partons cross the final state 
cut. 



3.2.1 One-parton contributions 

There are two different classes of contributions which have a single gluon crossing the final 
state cut, the real- virtual contributions and the contributions proportional to the LO r w 
soft function, derived by expanding the charge renormalization constant to 0(a s ). In what 



follows, Sft'^l 



, fi) and S R nt ' (A, fi) are simply the 0{a s ) terms in Eqs. ( |3.8| ) and 



respectively. 

The real- virtual interference terms can be derived by judiciously combining the analo- 
gous result derived in Ref. [24| for the hemisphere soft function and the all-orders- in-e LO 



results collected in Section 3.1. We find that the result can be written as 



S 



R-V, 
R I 



,k R ,\,n) = S% A (k L ,k R )r 2 * 5(X) + C A S R ut > {1) (X,fi) 

7rr(2 + e)r(l - 2e) cot(7re)r(-e) 2 
X r(l-e)r(-2e)(l + e) ' 



S(k L )S(k 



R 



>2r 



(3.11) 



where SX (&£, k R ) is defined in Eqs. (20) and (21) of Ref. [24]. It is worth pointing out 
that the in-jet terms, Sq (ki,k R ) r 2e (5(A), are simply derived by appropriately rescaling 
the real-virtual hemisphere integrals. For the sake of completeness, we recapitulate the 
result for (kL,k R ) derived in Ref. [24]: 



Sc A (k L , k R ) = C F C A fi 




4 2?H 32( 3 



7T 

: 30 



+ o( e 2 )) . 



(3.12) 



Given the results of Section 3.1, it is trivial to write down the charge renormalization 



contributions to the NLO t u , soft function: 



^ cn (fc L ,%,A,^) 



e 



S 



in,(l), 
R 



(k L ,^S(k R ) + S^ {1) (k R ,^5(k L ))5(X) 



+S R ut ' (1 \x, f ,)5(k L )5(k 



R 



(3.13) 



where /?o 



l -kc A 



fTp is the leading order /3-function. 



3.2.2 Two-parton contributions 

We now turn to the contributions with two partons crossing the final state cut. The sum 
of all such contributions to the Tu soft function has the form: 
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where ^Ic A and ^n* are the integrands for the CfCa and C^n/Tp color structures, for 
which explicit formulas can be found in Eqs. (16) and (27) of Ref. f24|]. The statistical 
factor of 1/2! in front of Ic A has its origin in the fact that, for the CfCa color structure, the 
two-parton contributions have two indistinguishable gluons in the final state. As mentioned 
above, the two-parton phase-space cuts are much more complicated than the one-parton 
ones: 

F 2 (k L , k R , A) = (-2iri) 2 8 {k 2 ) G (k°) 5 {q 2 ) 6 (q°) 

Q q n e k J(k R -k + - q+)5{k L )5(\) + @l®*5(k L - k~ - q-)5(k R )6(X) 

+e k n e q n 5(k R - k+)5{k L - q -)5{\) + @ q n e*6(k L - k~)s(k R - q +)6(\) 

+Q k n Ql ut 5(k R - k + )5 (x - ^y^) 5(k L ) + Ql& k ont 5(kn. - (x - ^y^") 

+e|G^ ut «5(fc i - k-)5 (\ - ^— p 1 ) 8(k R ) + ele k out 5(k L - q -)s (a - ^y^ 1 ) S(k R ) 



,,k n q A /\ q + + k + +q +k~ 



(3.15) 



There are four independent cases to consider: 



1. The same-side in-in contributions, where both partons are in the same jet (the second 
line of Eq. (gig)). 

2. The opposite-side in-in contributions, where one parton is in the n jet and the other 
is in the n jet (the third line of Eq. ( |3.15| )). 

3. The in-out contributions, where one parton is in a jet and the other is outside of both 
jets (the fourth and fifth lines of Eq. ( |3.15| )). 

4. The out-out contribution, where both partons are outside of all jets (the last line of 
Eq. (jug)). 

We consider each class of contributions in turn. 

Naively, one might expect the calculation of the same-side in-in contributions to be 
challenging since the evaluation of the analogous hemisphere integrals was by far the most 



technically demanding part of the calculation described in Ref. [24]. Fortunately, it turns 
out that we can recycle the same-side in-in contributions to the hemisphere soft function 
and obtain the desired result for free. The result of interest can be derived by appropriately 
rescaling the analogous hemisphere integrals. For arbitrary R we find 



Sg (k L , k R , A, n) = ( WgL + ^ j M *r fc f (A) I C F C A 



4 22 /134 
^ + 3? + I -9- 
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27 
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+ 



9 J 



8 40 152 4vr 2 /952 20vr 2 176<C 3 \ m , 



, (3.16) 



where we have made use of Eqs. (18) and (29) in Ref. [24j. 

For the opposite-side in-in contributions, the result can be expressed as 



(fcX) 1+2e 6{x) 



(3.17) 



where j c A { z i r ) an d /„. (z,r) are functions with e expansions that begin at O (e°) . In 
other words, 

fc A (*, r) = fcl (z, r) + /g (z, r)e + (e 2 ) 

f nf (*, r) = f { "] (z, r) + (z, r)e + (e 2 ) . (3.18) 

It is worth pointing out that, in the limit R — > 0, these functions are suppressed relative 
to the other contributions discussed in this section. In order to obtain exact results for the 
scale-dependent contributions to the integrated t u distribution, we need several moments 
of the functions defined implicitly above: 

•(0) 



/g'(0,r) = 16Li a (r 2 ) 



/i>,r) = 0. 



f bl (M = -16Li 3 (1 - r 2 ) + 16Li 3 



+ 32Li 2 (r 2 ) In 



l_ r i 



.° in 3 (1 - r 2 ) - 161n(r) In 2 (l - r 2 ) + V In (l - r 2 ) + 16C 3 
3 3 



/i>,r)=0. 



1 /9f (0) f <- 
1 dz ( Z J C A \ 2-z 

^ 



2 _ : -/g (0,0 



-8Li 2 (r 2 ) ln(r) + 161n(2)Li 2 (r 2 ) + 



88Li 2 (r 2 ) . „, 
+ 8Li 3 r 2 



16r 2 ln(r) 



and 



1 dz / 2 /»y (2= 



2-z 



3(r 2 -l) 3(r 2 -l) 2 3 
_ 32Li 2 (r 2 ) 16r 2 



ln(r) In (l - r 2 ) 



+ 



3 3(r 2 -l) 

32r 2 ln(r) i 32 
3(r 2 -l) 5 



+ — ln(r) In (l - r 2 ) .(3.19) 
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The calculation of the in-out contributions proceeds similarly. However, in this case, 
the results are significantly more complicated. As we shall see, these contributions lead to 
terms in the integrated t u distribution that depend on ^2 m a highly non-trivially way. 
The result can be written as 

+ C F n f T F g nf (^-,r^j + A*) J, (3.20) 

where g CjK (z,r) and g nj . (z,r) are functions with e expansions that begin at O (e°). In 
other words, 

9c A (*. r) = gP A (z, r) + pg> (z, r)e + (e 2 ) 

9n f (*> r) = gf f (z, r) + gW (z, r)e + (e 2 ) . (3.21) 

In order to obtain exact results for the scale-dependent contributions to the integrated t w 
distribution, we need several moments of the functions defined implicitly above: 

167T 2 

9c\ (0. *) = -f" " 64Li 2 (r) - 64Li 2 (-r) , 
< } (0,0 = 0, 

#g (0,r) = -16Li 3 (r 2 ) + 128Li 3 (l - r) + 128Li 3 f^y) " 128Li 2 (-r) ln(r + 1) 

64 

+64Li 2 (-r) ln(r) + 128Li 2 (r) ln(l - r) ln 3 (r + 1) + 641n 2 (l - r) ln(r) 

3 



-y7T 2 ln(l-r)-96C 3 



< } (0,r) =0, 



jr 1 ^ (^a o + r) _ 2g p A (o, r) ) = _ 16Li3 P) 

-64Li 3 (r) + 64Li 3 (y^jj + 64Li 3 (y^jj + 64Li 2 (-r) ln(r) + 64Li 2 (r) ln(r) 

16(r 2 + l) 64r 2 ln(r) 64 , „. , , . . 352 , . . , . 

~ 3(r 2_ 1) + 3(r2 _ 1)2 " y + 1) + 32 ln(r) ln 2 (r + 1) + — ln(r) ln(r + 1) 

+ f vr 2 ln(r + 1) + ^ ln(l - r) ln(r) - 32C 3 - ^ , 

jf £ « CD + < <i„ r , - 2<) (M ) = -« - ^ 
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and 



dz 2 

Z2 Jo 



1 dzi 

Zl 



£2 



r J " < } (0, r) j = X n/ (x, r) - Xnf (1, r) . (3.22) 



The functions (x, r) and x n ^, (x, r) are given in Appendix [B] in terms of an appropriate 
set of one- and two-dimensional harmonic poly logarithms. 3 
Finally, for the out-out contribution, the result is 



S R (k L ,k R ,X,fi) 
176 ln(r 



4c 



(2A) 



l+4e 



5{k L )5{k R ){C F C A 



321n(r 



+ 321n 2 (r) + 64Li 2 (-r) - 64Li 2 (r) + 64Li 3 
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128Li 3 




128Li 2 (-r) ln(l 



r + 1/ \r + 1 

+ 128Li 2 (r) ln(l - r) - 64Li 2 (-r) ln(r) + 128Li 2 (-r) ln(r + 1) + 64Li 2 (r) ln(r) 



-128Li 2 (r) ln(r + 1) 



16 (r 2 + 1) 64r 2 ln(r) 



32 



3-3r 2 



3(r 2 
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ln 3 M 



641n 3 (r) 



— ln 3 (r + 1) + 1281n(r) ln 2 (l - r) - 32 ln(r + 1) ln 2 (l 



32 
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176 ln 2 (r) 
3 

1072 m(r) 
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-CFUfTp 



+ 641n(r)ln 2 (r + l) 
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ln(r) ln(l 



)-321n 2 (r + l)ln(l-r) 
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7r 2 ln(l - r) H 7r 2 ln(r) 
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— ln(r) ln(r + 1) + — ir 2 ln(r + 1) 



641n(r) + 256Li 2 (r) + 32(1 + r 2 ) + 128r 2 ln(r 



9 

641n 2 (r) 



-— — ln(r) In 



3c 
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r + 1 



256 W X , / 

H m(r j m(r + 1 

3 



3(l-r 2 ) 3(1 
128vr 2 



9 



-r 2 ) 2 



+ 



3201n(r) 
9 



(3.23) 



The calculation of this contribution is actually quite involved; in the integrals that lead to 
Eq. ( |3.23| ), there is explicit R dependence (which does not factorize from the rest of the 
expression) and one has to use sector decomposition 5C] in a non-trivial way to deal with 
the phase-space singularities. As we shall see later on, integrating the out-out contribution 
to the Tu soft function leads to logarithms of We expanded the non-trivial part of the 
right-hand side of Eq. ( 3.23j ) to O (e°) because it turns out that this completely determines 
the logarithms of which appear in the integrated distribution; the general arguments 
of Q forbid terms that go like In 3 (£) or In 4 (£). 

Let us now briefly summarize what we have accomplished. In this section, various 
contributions to the NLO t u soft function were calculated. The full result, S R (k^, k R , A, //), 
is naturally written as a sum of six contributions: 



S$\k L , k R , A, y) = S^~ v (k L , k R , A, y) + S% cn (k L , k R , A, M ) + S%(k L , k R , A, n) 



?R-V 



Ron / 



3 For the reader less familiar with harmonic polylogarithms, we recommend reading Ref. 
pendix A of Ref. G§. 



and Ap- 
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+S R 2 (k L , k R , A, fj.) + S T i(k L , k R , A, p) + S%(k L , k R , A, fj). (3.24) 



The three terms in the first line of Eq. (3.24) were easily calculated using known results 



(computed in Refs. and J24j). The three terms in the last line were more involved. 
Instead of computing S T ^(kL,k R , A, fi) and (k^, k R , A, fi) directly, we instead computed 
the moments of those contributions relevant to the integrated jet mass distribution. We 
also chose not to compute the O (e) coefficient of S T ^(kL, k R , A, fi), which only affects the 
constant part of the integrated distribution. The results given in this section determine 
the unknown scale-dependent terms in the integrated r w distribution in the threshold limit 
at NLO. In the next section, we examine and check these results. 

4. The Jet Thrust Distribution at Two-Loops 

In the previous section, the parts of the two-loop t w soft function relevant for the differen- 
tial jet thrust distribution were calculated in dimensional regularization. We have chosen 
not to compute the parts of the soft function that contribute only to the zero bin (r w = 0) 
of the differential jet thrust distribution. For the hemisphere soft function, these terms 
were calculated (see Ref. j^4|) because they are important for precision a s measurements. 
However, in this paper, we are primarily interested in understanding the structure of the 
non-global contributions to the integrated distribution and in testing the refactoriza- 



tion conjecture of Ref. [30]. We therefore followed the approach of Ref. [33| in Section 
|3| and calculated only those contributions to the soft function which depend on ratios of 
dimensionful scales. 

The triply differential jet mass cross section is a complicated singular distribution which 
does not have a simple expression in momentum space. Although the two-loop integrated 
jet thrust distribution is not simple either for finite R, it is quite a bit easier to work with 
in practice due to the fact it is free of distributions. In terms of the various contributions to 
the Tu soft function studied in Section ||, we can write the integrated jet mass distribution 
as 

Kr\t w ,u),h) = J^ dr' u dA dk L dk R S^\k L , k R , A, fj)5 (r' u _ fcx + 

= / dr^ / dA / dk L dk R S^- v (k L ,k R ,X,fi) 
Jo Jo Jo L 

+S% en (k L , k R , A, /x) + (k L , k R , A, p) + S%(k L , k R , A, fj) 
+ S^(k L ,k R A,p) + S^(k L ,k R A,p)]5^-^±^^j , (4.1) 

where the superscript (2) reminds us that we are calculating the integrated distribution 
at NLO. The scale-dependent parts of the six terms in this expression are given explicitly 
in Appendix When these soft contributions to the integrated r w distribution are ap- 
propriately combined with the contributions coming from the hard and jet functions, the 
result is a concrete prediction for the integrated distribution in the threshold region. 
This prediction should be valid for arbitrary R, up to power corrections in t u and ui/Q. 
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Figure 2: The difference between the SCET and the QCD caiculations of ^37- for various values 
of R and uj. The left plots show the CfCa color structure and the right plots shows the CpUfTp 
color structure. The top plots show the variation with R at fixed ui = 10 _2 Q and the bottom plots 
the variation with uj at fixed R = 1CP 2 . For the QCD calculation, a cone algorithm is used with 
the jet direction taken to be the thrust axis. 



4.1 Comparison to full QCD 

We would like to compare our analytic results for the threshold limit to full QCD, for which 

(2) 

the distribution is only known numerically. To do so, we first differentiate K R (t u ,uj, /i) 
with respect to t w to determine the NLO part of the differential distribution. This result 
is then appropriately combined with the LO contributions coming from the hard and jet 
functions, as prescribed by the factorization formula, Eq. ( |2.2| ). The resulting expression 
is the threshold prediction for y--^-- This prediction can then be expanded to O (a^) 
and numerically compared to QCD using the program event2. The details of the a s 
expansion can be found in e.g. Ref. 

The comparison shown in Figure [2] is the difference between the full QCD distribution 
calculated using event2 and the two-loop prediction using SCET, with for various values 
of R and uj. The plots show that SCET is reproducing all of the singular dependence in r w , 
for all values of R up to power corrections in uj/Q. This is a non-trivial check that SCET, 
through factorization formula (p.2j), is reproducing the correct singular behavior in r w for 
any R and sufficiently small uj. Note that one needs t w -C tj to approach the threshold, 
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so for very small u, the threshold migrates into a numerically unstable regime, which is 
apparent in the figure. 



4.2 Analytic checks and definition of the non-global contribution 

(2) 

Next, we can compare the //-dependent parts of K R (t w ,u;,/x) to the \i dependence pre- 
dicted by the factorization theorem in SCET (Eq. (j2.2j)). With the one- and two- loop 
anomalous dimensions for the jet and hard functions, the one-loop QCD /3-function coef- 
ficient, and K^ R (t^ , uj , //) as inputs, SCET predicts the \i dependence of the integrated t u 

(2) 

distribution exactly. As it must, the \i dependence of K R (r w ,w, /u) exactly matches the 
prediction of the factorization theorem. 

In order to say more about the non-global contributions to the integrated jet thrust 

(2) 

distribution, it is convenient to subtract off the //-dependent part of K R (r w ,w,/u). In fact, 
as we will argue in Section || the preferred way of doing this employs the integrated form 
of the refactorization ansatz reviewed in Section ||. In integrated form, the refactorization 
formula of Eq. (2.8) is 

K R (T u ,u,fi) = i^ n ( W )i^ ut ( W , A) ® K f R {r u ,U3) . (4.2) 

In order to make use of the above formula, one also needs the anomalous dimensions 
associated to i^j^r^,//) and K R nt (uj,fi), 

J$ = 75 + T cusp In (j^x) and ^ = " r - s P hl (iTr) ' (43) 



In Eqs. ( |4.3|) above, 75 is the thrust soft function anomalous dimension. Eq. ([4.2;), 
Eqs. ( |4.3| ), and some of the ingredients mentioned in the previous paragraph determine a 
function ^ efac (r w ,w , fi) which one can use to define the non-global contributions to the 
integrated t u distribution. Explicitly, we have 



= CfCa 



two— loop 



{jL) /-_176Li 2 (-r) _ 44 _ 8 5361n(r) _ 44^ 
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The difference of the exact result and (t^, u, h) is 
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176Li 2 (r 2 ) . „. 
4^ + 8Li 3 (r 2 ) 
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64r 2 ln(r) 
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+ ^ln(r)ln (l - r 2 ) ln 
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3 9 / V 2a> 



3(l-r 2 ) 3 (l-r 2 ) 2 

(4.5) 



As required, this difference is //-independent and contains all of the non-global logarithms. 
The functions Xc A ( x > r ^ anc ^ X nf ( x ' r ) m Eq. (L5) are defined in Appendix [B|. There 



should be an additional t w - and w-independent function of R in Eq. (4.5), which has not 
been calculated. As mentioned above, this function would not affect the differential t w 
distribution and is therefore of secondary importance. 

In recent years, quite a bit of effort has been devoted to understanding the structure of 
NGLs. So far, most of the analytical calculations have focused on the leading logarithms. 
For example, in Eq. ( |4.5| ) above, the leading NGL is simply 



O r i(l6Li 2 (r 2 )-^ In 2 



(4.6) 



which agrees with results from Refs. [30, |33j. Leading NGLs are certainly of interest but 
such computations are not, by themselves, likely to be of much use if our ultimate goal 
is to understand the resummation of NGLs. On the other hand, our exact result for the 
non-global contribution to the integrated t u distribution is extremely complicated (both 
Xc A ( x ' r ) anc ^ Xn f ( x > r ) depend on two independent scales in a highly non-trivial way) 
and this makes it difficult to gain any further theoretical insight. If there is any hidden 
structure in our results, it should be significantly easier to identify in the small R limit. 



5. The Small R Limit of the Non-Global Logarithms 

It was observed in Ref. that, in the small R limit, the jet thrust soft function simpli- 
fies. Having computed the exact scale-dependent contributions to the integrated jet thrust 
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Figure 3: By rescaling one of the jet directions a small cone can be mapped to a hemisphere. 
Rescaling in this way makes some universal properties of non-global logarithms easier to identify. 



distribution at finite R, we can now understand more precisely what is happening. As we 
will see, in the small R limit, the non-global contribution to the integrated t u distribution 
is determined completely by the non-global contribution to the integrated hemisphere soft 
function. 

It has been noted by a number of authors that, for jet observables defined with the 
anti-fcr jet algorithm, one typically finds that the R — > limit of the leading NGL is simply 
given by the result in the hemisphere case, up to possibly a factor of two depending on 



the observable in question [29, |3C], ^1], 33]. The jet algorithm employed in this paper is 
very closely related to the anti-/c^ algorithm. In fact, the in-in and in-out contributions 
are exactly the same in both jet algorithms. The only difference between the algorithms 
at two loops is that they treat the out-of-jet radiation differently. We strongly suspect 
that the NGLs in the two jet algorithms actually coincide, at least in the small R limit. 
This implies that we should be able to check whether, for the integrated anti-fcr jet thrust 
distribution, one can relate the next-to-leading NGLs to analogous NGLs in the integrated 
hemisphere soft function calculated in Refs. |3(]] and since we have computed the exact 
expression, we are now in a position to go beyond the leading-logarithm approximation. 
In what follows, we assume some familiarity with the conventions and definitions that we 
established in Section ||. 

5.1 In-in rescaling 

Before discussing the in-out contributions and non-global logarithms, it is instructive to 
first consider the simpler case where both soft particles get clustered into the same jet. 
The same-side in-in contribution to the integrated distribution is 

2 e^\ 2t f d d q d d k 



(2ir) d {27r) d 



x I ^Io A + InA (-27ri) 2 6(q 2 )Q [q°) 5(k 2 )Q (k°) 9(rAr - k+) 
xQ(rq- - q + )5(k L )5(k R - q+ - k + )5 (A) 5 (i> u - , (5-1) 

where Iq a and I nf are the two-loop, two-parton integrands for the CfCa and CprifTp 



color structures, given respectively by Eqs. (16) and (27) of Ref. 24]. 
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To simplify the theta and delta functions, we can rescale the n-jet Wilson line. This 
turns the cone jet into a hemisphere jet, as illustrated in Figure 5.1. At the level of the 
integrand, this rescaling amounts to making the replacements: 

k~ — > and q~ — > . (5.2) 

This change of variables induces the transformations 

Ic A -> r 2 Ic A and I nj -)• r 2 I nf (5.3) 

and the phase-space cuts in Eq. ( |5.1| ) that depend on the minus components of the light- 
cone momenta become 

0(AT - k + )@(q~ - q + ). (5.4) 

Remarkably, the transformed integrand reduces to the O (a 2 .) same-side hemisphere inte- 
grand times an overall factor of r 2e . A similar rescaling was observed in the context of 
beam functions in [51]. 

After performing the trivial integrations over ki and k^, the result is simply expressed 
in terms of the same-side contributions to the integrated hemisphere soft function: 



K%(t u ,u>,h) = 2r 2e H QAr'^ f" dA 
Jo Jo 



S(X)n A ' 

l+4e 



{QtL 



C F C A g h ST (<0 + C F n f T F g^f (e) 



(5.5) 



where Qq^ 11 ^) and <7^ mi (e) are the same-side contributions to the NLO hemisphere soft 
function, given in Eqs. (18) and (29) of Ref. Q. 

We used this rescaling to compute the two-loop same-side in-in integrals in Section ||. In 
fact, this result is straightforward to generalize. At any order in perturbation theory, if all 
of the particles go into a single jet, the contribution to the integrated t w soft function will be 
given by the analogous hemisphere calculation with the light-cone momenta appropriately 
rescaled in the manner described above. 

Before moving on, we briefly consider the case where the two partons end up in separate 
jets and show that this contribution is suppressed in the small R limit. First, for any r, 

(2) 

this contribution to K R (t u ,uj, y) is given by 

T " , / r f°° „ „ /„ ^4/Ve 7E V e f d d q d d k 



K%(t u ,oj,») = 2 £ dri £ dA £ dfc L dM4vr) 4 (^) / 



(2?r) d {2ttY 

1 



2 -I Ca + l nf {-2m) 2 5{q 2 )Q [q°) 5(k 2 )Q (k°) Q(rk~ - k+) 
xQ(rq+ - q~)5{k R - k + )5{k L - q~)S (A) 6 (r' u - • (5.6) 

Let us try to proceed as above, this time rescaling k~ and q + : 

k~ , q + 
kT — > and q + — > . (5-7) 
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This time Ic A and I nf do not transform simply under the rescaling. Instead we find that 

Ic A ^r 3 I CA {r) and I nj -> r 5 I nf (r) , (5.8) 

where Ic A ( r ) an d In s { r ) are regular but non- vanishing at r = 0. Since, as before, the 
integration measure contributes a factor r _2+2e , we see that this entire contribution is 
suppressed at small R relative to the one studied above. In fact, in the small R limit, the 
opposite-side in-in contributions are suppressed relative to all of the other contributions 
and can therefore be neglected. 

5.2 In-out rescaling at small R 

Next, we will look at the contributions where one soft parton gets clustered into a jet and 
the other does not. Here we will find a similar mapping to the hemisphere integrals but 
only at small R. At the level of the integrand, the in-out contribution to the integrated 
distribution is 



-^R 3 (-r w ,w,/i) 
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(5.9) 



Changing variables in an attempt to map the cone jet to a hemisphere jet as above, 
we see that the phase-space cuts in Eq. ( |5.9| ) that depend on the minus components of the 
light-cone momenta can be written as 

2r@(k- - k + )Q(q- - r 2 q + )Q(q + - q-)5 (q + r + q~ - 2r\) . (5.10) 

This is not exactly a hemisphere projection for general R. However, in the small R regime, 
the above expression simplifies to 

2r9(fc- - k + )9{q-)0(q + - q~)5 (q~ - 2r\) . (5.11) 

The phase-space cuts enforced by ( |5.11| ) are identical to the ones which arise in the 
calculation of the opposite-side contributions to the two-loop hemisphere soft function, 
provided that one makes the replacement ki — > 2r\. To see this, after performing the 
trivial integrations over kz and k R in Eq. (|5.9|) above, replace the phase-spaces cuts that 
depend on k~ and q~ using ( |5.1lD and then compare the resulting product of thetas and 
deltas in the integrand to the second line of Eq. (13) in Ref. [24| (with ki = 2r\ and 
= Q T D- With this understanding, we find that the above expression can be written in 
terms of the opposite-side contributions to the integrated hemisphere soft function: 
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(5.12) 
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q a (z, e) and J (z, e) are the opposite-side contributions to the NLO hemisphere 
soft function, defined in Eqs. (17) and (28) and Appendix A of Ref. [24]. In deriving Eq. 



(5.12) we found it useful to make the change of variables 



x 

Q 



A 



y_ 

2r 



(5.13) 



This correspondence between the in-out contributions to the integrated distribution 
and the opposite-side contributions to the integrated hemisphere soft function is valid for 
sufficiently small R and arbitrary values of and w. The striking similarity between Eqs. 
Q5.5| ) and (5.12) is suggestive. In fact, for sufficiently small R, we will show in the next 
section that it is possible to very accurately model all of the terms in the integrated r w 
distribution not captured by the refactorization ansatz - including all of the NGLs - using 
the non-global contributions to the integrated hemisphere soft function. 

5.3 Small R limit 

We have shown that the in-out contributions at small R (and so r ~ R) are given by 
the opposite-side contributions to the integrated hemisphere soft function times an overall 
factor of 2R 2e . We were inspired by this small R correspondence to take a closer look at 
the non-global contributions to the integrated hemisphere soft function, IZf (z). 

we are able to take the small R limit of the 



Using the results of Sections 5.1 and [57 
scale-dependent non-global contribution to the integrated t w distribution that we presented 
in Eq.(p|). We find 



A' 
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(r w Q\ 
\2Rjjjj 



(5.14) 



defined in Eq. (53) of Ref. [24]: 



where TZf(z) is the non-global contribution to the integrated hemisphere soft function, 
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(5.15) 



9 27 

In deriving this small R limit, the rescaling of the in-out integrals in particular was critical 



- it is certainly not straightforward to derive Eq. ( 5.14j ) directly from the expression for 
A^(r w , uj) given in Eq. (|4.5| ). Again, there may be additional contributions to K^^^t^jU) 
which are r w - and ^-independent and have not been computed. 
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Figure 4: Contributions to the differential jet thrust distribution coming from the non-global terms 
in the soft function. We show ;nn7rT '^fl(' r ) ( souc and its small R limit, j^if^Jr) (dashed), 
plotted as a function of r = for various R. 

In Figure HI we compare the r w distributions coming from the exact expression, kL(t w ,uj), 

£ 

and from the small R approximation, K j r .q{t u ^u), for both non-trivial color structures. 
We see that the approximation works spectacularly well and remains valid for small but 
experimentally relevant values of the jet cone size (R ~ 0.1). 



5.4 Extraction of the non-global logarithms 



£ 

Traditionally, one would define the non-global logarithms as those terms in K R (Tu, to) which 



diverge in the 



In 



2uj 



3> 1 limit. However, the non-global logarithms in our integrated 



jet mass distribution are somewhat more subtle since there is a term in K^(t u ,uj) which 
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diverges in the R 



limit as well. To sensibly extract the NGLs in our case, we actually 

We found in 



(5.16) 



need to consider the double limits 1 > R > and 3gf > 1 > R 
Section [5l| that 

JSr£(r ttl u) = 2% (gg) + /(H) + O(iJ) , 

where f(R) is a w-independent and r^-independent function of R which we have not cal- 
culated. The appearance of the function TZf(z) in Eq. ( |5.16| ) suggests that the leading and 
next-to-leading NGLs in the integrated t u distribution ought to be simply related to the 
leading and next-to- leading hemisphere NGLs derived in Refs. [24, plf . 



In 



We can easily extract the NGLs in the integrated distribution by first taking the 
(it?) ^ 1 limit of the non-global function K^ir^^oo) and then taking the small R limit 

of the expression that results. In this approximation, we find that the function K^(t w ,o;) 
reduces to 
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Strictly speaking, Eq. ( 5.17 ) is only correct up to terms independent of r w and uj (i.e. 
functions unconstrained by the calculations performed in this paper). However, although it 



remains to be conclusively proven, we strongly suspect that Eq. ( 5.17 ) actually captures all 
of the singular R dependence present in the contributions to the integrated r w distribution 
not determined by the refactorization ansatz. After multiplying the above expression by 1/2 
and making the substitution In (jjrgj^j — > hi(z), we find that Eq. ( 5-17] ) exactly reproduces 
the hemisphere NGLs of Eq. (54) in Ref. |2^1 . This shows that, as hoped, the next-to- 
leading NGLs in our integrated jet thrust distribution are simply related to corresponding 
next-to-leading NGLs in the integrated hemisphere soft function. In fact, if it turns out 
that, as seems likely, the anti-/cr NGLs coincide with ours in the small R limit, we have 
shown that there is no essential difference between the small R NGLs that arise in the 
integrated anti-Zc^ jet thrust distribution and the NGLs that arise in the hemisphere mass 
distribution. 

Clearly, the simplicity of the NGLs derived above suggests that a simple factorized 
picture emerges in the small R limit. It seems likely that, for sufficiently small R, the 
soft function refactorizes to all orders in perturbation theory: 



SR(k L ,k R ,X,fi) = Sn 



RH 



Si, 



k R \ ( 2VRX n 

'-'out I > n 



y/R^L 



(5.18) 



Here, the functions S^z) and Sf(z) denote respectively the ^-dependent and /i- independent 
parts of the hemisphere soft function |l^, |24], |37j . Writing the soft function in this refac- 
torized form makes the relationship between the small R exclusive jet mass soft function 
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and the hemisphere soft function completely explicit. It is worth pointing out that S^(z) 
and Sf(z) have no explicit R dependence; their R dependence is simply encoded in their 
arguments, -^7=^ and in fact, the only part of ( 5.18| ) which is not a simple rescaling of 



part of the hemisphere soft function is 5 ou t f 2v ^ A ; Rj ■ Moreover, the dependence on the 
first argument of S out (^^,It) is fixed by RG invariance. Actually, it is possible to fix 



the entire singular R dependence of S out (z, R) with arguments similar to the ones made in 



Sections |5.l| and [^2} This analysis, however, is beyond the scope of the current paper. 



When we expand Eq. (|5.18[) to O (a 2 s ), we find that 



K' r (t u ,u,) = 2K f + O(R) . (5.19) 



This is completely consistent with the function K^_ >0 (r^, cj), Eq. ( |5.14 ), coming from the 



expansion of the exact result. It also predicts that the unknown r w - and w-independent 
function of R in the integrated distribution should vanish in the small R limit. Furthermore, 
since we have confirmed this refactorization at two loops, including non- global logs, the 



refactorization and r cusp ansatz of Ref. |3(| are also confirmed at two loops - Eq. ( 5.1£j ) 
reproduces them when Sf(z) = 1. 

6. Conclusions 

Over the last several years it has become increasingly apparent that jet substructure will 
play an important role in collider physics, and that theoretical calculations of jet sub- 
structure will require new organizational tools. For example, the simplest substructure 
observable, jet mass, when examined in an exclusive context, such as for the hardest or 
two hardest jets, is associated with poorly understood dependence on the multiple relevant 
scales (jet masses and veto scales). In this paper, by studying in detail the non-global 
structure of the soft, scale-dependent contributions to the integrated cone jet thrust dis- 
tribution at two- loops, we took an important step towards understanding the interplay of 
these scales. We presented compelling evidence that the non-global logarithms that arise 
in a realistic class of jet algorithms when expanded around small jet size R are in fact no 
more complicated than the non-global logarithms that arise when one uses the hemisphere 
jet algorithm. 

It seems that at each order in perturbation theory, it is possible to take the limit of small 
jet cone size in a consistent and useful way at the level of the integrands and that taking 
this limit effectively reduces all subsequent integrations to the corresponding hemisphere 
ones. This vastly simplifies the non-global parts of the distribution. A comparison of the 
non-global contribution to the integrated hemisphere mass distribution given in Section 
|5.3| and the analogous integrated jet thrust distribution given in illustrates the point - 
the hemisphere results of Ref. [^J] reproduced in Eq. ( |5.15 ) are far simpler than the finite 



R results in Eq. ( |4.5[) , derived in this paper. Naively, one might think that this procedure, 



explained in more detail in Sections 5.1 and 5.2, only gives an accurate approximation if 



one works with unrealistically small jet resolution parameters. However, as illustrated in 
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Figure ||, this is not the case at all; we found that approximating our result for the finite 
R differential distribution using this procedure works great even for R ~ 0.1. 

One point worth emphasizing is that the refactorization ansatz of Ref. p0[j played an 
important role in our analysis. In Figure [l] we compared the O (a^) prediction of this 
ansatz to full QCD at the same order using event2 and showed that the refactorization 
ansatz reproduces the singular terms in the QCD differential distribution very nicely, up to 
subleading jet algorithm and NGL effects. Indeed, in Section ||, we defined the non-global 
contribution to the integrated jet thrust distribution by subtracting the refactorization 
ansatz in integrated form from our complete result for the scale-dependent part of the 
integrated distribution. If we had used a different subtraction procedure it is not clear that 
we would have been successful in making contact with the integrated hemisphere function. 
This remark is related to the fact that, at the outset of this work, it was not clear what 
appearance, if any, the jet cone size R would make in the non-global contribution that 
remains once one subtracts off the integrated refactorization ansatz. 

Our results suggest that the non-global contribution to the integrated jet thrust dis- 
tribution at small R is simply 



xkT*iT u rf = 2K f (^ , (6.1) 



where TZf(z) is the non-global contribution to the integrated hemisphere mass distribution. 
We have only shown this up to a t u - and w-independent function of R, due to the fact that 
we consistently dropped all contributions that depended on no scales at all or on R but 
not on a ratio of dimensionful scales. This small R limit depends on the precise definition 
of K^TuiiUj). In our definition, K* r (tu : u)) is everything not given by the refactorization 
from [ [30f| . With this definition, Eq. ( |6.lD implies that all of the remaining singular R 
dependence is inextricably linked to the ratio ^gfi. 



The assertion of the last paragraph is very sensible because we showed in Section 5.4 
that the leading and subleading non-global logarithms are 



kST^> - c F c A _ £ * (g) + ( - 1 + ^ - 16(3 ) y (£) 

32vr 




\2Roj J 



+ ■■■. (6.2) 
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In the appropriate limit, Eq. ( |6.1| ) reduces to Eq. ( |6.2[ ) up to terms not constrained 
by our calculation. The form of Eq. ( |6.2j ) also shows that the NGLs that arise in the 
integrated cone jet thrust distribution are naturally R dependent. The great strength of 
our approach is that we abandoned the leading-logarithm approximation, preferring instead 
to first perform exact calculations and only then attempt to determine an appropriate 
approximation scheme. 

The results recapitulated above motivated us to modify the refactorization formula of 



Ref. [Bel] to include the non-global contributions as well. In Section 5.4, we discussed the 
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possibility of a refined refactorization formula 



Sn(k L ,k R , A,/t) = S/j, ( -%-] f-^=-\ S out ( 2 ^^ ,R 



9 ^(^) s '(^j + ^' (6 - 3) 

where the functions S^z) and <S/(z) denote respectively the //-dependent and //-independent 
parts of the hemisphere soft function |T^, j37| . Eq. ( |6. 3[) is consistent with all of the obser- 
vations made in this paper and, if true, would make precise the statement that resumming 
the NGLs in the cone jet thrust distribution is equivalent to resumming the NGLs in the 
hemisphere mass distribution. 

In future work, it would be useful to understand to what extent refactorizations like Eq. 
Q6.3D hold beyond two loops. If this is understood, similar arguments should be appropriate 
for precision predictions of jet substructure observables at hadron colliders. In any case, 
it is of great interest to try and generalize Eq. ( |6.3| ), both to multi-jet processes at e + e~ 
colliders and, if possible, to di- and multi-jet processes at hadron colliders. 

In order to understand non-global structure in more detail, it may be necessary to 
explore jet algorithms distinct from the thrust cone algorithm used in this paper. Most 
prominent among the alternatives is the ant i- fey jet algorithm which, as noted in Section 
H is very similar to the thrust cone one and has the advantage that it is widely used by 
experimentalists. It would also be interesting to have all of the i?-dependent constant 
terms that were dropped in writing down the non-global contribution to the integrated jet 
thrust distribution, since they would significantly clarify the structure of the ln(i?) terms 
that remain in the small R limit. 

Finally, even if Eq. (|6.3|) turns out to be a useful approximation, one still has to deal 



with the hemisphere NGLs themselves. Resummation of these NGLs seems difficult pO, 22 



24, 25, 2qj . Perhaps one can choose scales so that the non-global structure has a numerically 
small effect. Nevertheless, it would be nice to understand NGLs in more detail, beyond 
the leading NGL. One way forward would be to perform an explicit three-loop calculation 
of the integrated hemisphere soft function. While it goes without saying that this will be 
no easy task, we expect that the integrals which arise at three loops will turn out to be 
solvable using modern multi-loop computational techniques. 
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A. Calculation of the Integrated Jet Thrust Distribution 



In this appendix, we calculate the six terms in Eq. ([O]) in the order in which we discussed 
the corresponding contributions to the NLO soft function in Section ||. We begin with 
the charge renormalization contributions. The result is 



K^(r m uj,n) = C F C A 
22tt 2 



176 
~9~ 



ln d 



^-)+^ln(r)ln 2 ( JL 



+ 



441n 2 (r) 



9 



,n | Jt-\ + l n (JL\ (^U 2 (-r) 441n^) 44^\ _ 88 

t w QJ \2uj)\ 3 3 9 J 3 KI \2u 



+CFUfT F 



64 i 3 

— ln d 



32 



— ln(r) In 2 , . , , 

tM) 3 V ; \tM) V 9 



A* 



^8vr 2 161n 2 (rj^ ln / // 



+ ln(— ) r- 64Li2 '- r ) - 15 ln» - ^ + » ln W In' f i 
V2w/ I 3 3 w 9 / 3 v ' \2ui 



(A.l) 



For the real- virtual interference contributions, the result is 
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and for the same-side in-in contribution, we have 
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The calculation of the opposite-side in-in contributions is somewhat less trivial since 
the integral over (ki, A, /x) requires sector decomposition. Carrying out this analysis 
leads to 
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for the finite part of the scale-dependent terms. Eq. ( [A. 4 ) is written in terms of the 
functions introduced in Section |3[ We can easily evaluate the above expression using the 
moments tabulated in Eqs. ( 3.1 9| ) and the result is 
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The calculation of the in-out contributions requires sector decomposition as well. In 
terms of the functions introduced in Section |||, the scale-dependent terms are given by 
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Using the moments tabulated in Eqs. ( |3,22 ), we find an explicit expression for Eq. ( [A. 6 ): 
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where Xc A r ) an< ^ Xrif ( x > r ) are non-trivial functions built out of one- and two-dimensional 
harmonic polylogarithms. They are given explicitly in Appendix [B]. 
Finally, for the out-out contribution, the result is 
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Now that all the pieces are in place, we can combine them together and study the result. 



B. Analytic Expressions For Non- Trivial Two-Parameter Integrals 

The non-trivial two-parameter integrals 

first encountered in Section ||, were parametrized there in terms of two functions, \ Ca {x, r) 
and X n ( x ' r )> such that: 

As explained in Section ^, C/^(x,r) and gffl(x, r) are simply the leading order terms in 
the epsilon expansion of the in-out contributions to the r w soft function. In this appendix 
we provide analytical expressions for both Xc A ( x,r ^ an< ^ Xn/( x > r )- ^ n w hat follows the H 
functions are one- or two-dimensional harmonic polylogarithms (introduced in Refs. |48[| 
and [jljj respectively). 

For the CpCa color structure we have 
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